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ABSTRACT 

We show that the candidate progenitor of the core-collapse SN 2011dh in M51 (8 Mpc 
away) was fading by 0.039 ±0.006 mag/year during the three years prior to the supernova, and 
that this level of variability is moderately unusual for other similar stars in M 51. While there 
are uncertainties about whether the true progenitor was a blue companion to this candidate, 
the result illustrates that there are no technical challenges to obtaining fairly high precision 
light curves of supernova progenitors using ground based observations of nearby (< 10 Mpc) 
galaxies with wide field cameras on 8m-class telescopes. While other sources of variability may 
dominate, it is even possible to reach into the range of evolution rates required by the quasi-static 
evolution of the stellar envelope. For M 81, where we have many more epochs and a slightly 
longer time baseline, our formal 3a sensitivity to slow changes is presently 3 millimag/year 
for a My ^ -8 mag star. In short, there is no observational barrier to determining whether the 
variability properties of stars in their last phases of evolution (post Carbon ignition) are different 
from earlier phases. 

Subject headings: supemovae:general, supemovae: individual: SN 201 Idh 



1. Introduction 



The last few years have seen stead y progress in the identification of the progenitors of core-collapse 
supernovae (ccSNe, see the review by ISmartti 120091) . The progenitors of Type IIP SNe are red super- 
giants, although there is some evidence that the most massive progenitors are less rnassive than the ex- 
pected upper mass of red supergiants at the ends of their lives (IKochanek et al.ll2008l ISmartt et al.l 12009 , 
but see IWalmswell & Eldridgel l201lh . The progenitors of t wo Type lib supemo vae have been identi- 
fied, a (pro bably) mass transfer binary system in S N 1993 J (|Aldering et all 1 19941) and the progenitor of 



SN 201 Idh (Maundetal 
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l|). In two cases, SN 2005gl (IGal-Yam & Leonardll2009h 



and SN 1961V ( Kochanek et al.ll2011l. iSmith et al.ll201lh . the progenitors appear to be massive, luminous 
blue variables (LB V). The dust-e nshrouded progenitors of Type Iln SN 2008S ( Prieto et al.ll2008l) and 2008 
NGC 300 transient (|Prietoll2008i) were identified in the mid-infrared. No progenitors of Type Ibc SNe have 
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been identified, and tiiis will likely remain challenging due to their lower rates, lar ge expected bolo metric 
corrections, and the ease with which a binary companion can be optically brighter (IKochanekl 120091) . Still, 
progenitor studies are now well established. 

The next frontier is the variability of progenitors. Very little is known observationally about the vari- 
ability of these stars shortly before explosion. The proge nitor of the Type Ilpec SN 1987 A varied by less 
than a few tenths of a magnitude during its last century (see lPlotkin & Claytonll2004l and references therein), 
while the progenitor of the Type lib SN 1993J varied by less that 0.2 mag over a 6 month period 9 years 
before the explosion dCohen et al.lll995h . The progenitor of the Type IIP SN 2008cn was probably vari- 
able at a level of ~ 0.2 mag and the sparse light curve is potentially interpretable as an eclipsing binary 



(lEhas-Rosa et al.l 120091). The pr ogenitor of the Type lb SN 2006jc showed an outburst two years before 
explosion /Pastore llo et al. 20071 ) . SN 196iy also had an outburst prior to explosion if it is interpreted 
as an SN (Kochaneketal.lboi;!, ISmith et al.ll2011h . In several tr ansients, variabilit y was either ob served 
(SN 2010da, .Laskar et al. 20 1Q) or well-constrained (SN 2008S. IPrieto et all boosl: NGC 300-OT, |prieto 
20081 . [Thompson et al.ll2009l) in a dusty stellar wind rather than directly in the star. More generally, out- 
bursts in the last ^centu ry before explo sion may be required to explain the post-explosion evolution of 
many Type Iln SNe (e.g.. lFox et al.ll201ll) . In short, data are available for very few SNe, and when it exists 



it is generally too sparse to interpret. 

Theoretically, all SN progenitors are variable because their envelopes are evolving on a thermal time 
scale in response to the rapid changes in the core luminosity. In the standard lore, this quasi-equilibrium 
evolution is too small to be observable. Rough estimates can be e xtracted from some t abulated evo lution 
models, and typical ra tes in the optical are 0.1-1.0 millimag/year ( Schaller et al. 1992, Heger et aPEoOQ . 
Heger & Langen 120001) . but none of these studies was really intended for studies of the surface evolution 
during the last century. While these rates are small, it is certainly possible to achieve the necessary photon 
counting statistics to detect 1 millimag/year changes in progenitors for nearby (< 10 Mpc) galaxies (see the 
discussion in ®. 

It is likely in most cases that other sources of variability will dominate and mask the slow evolution 
of the envelope. Many massive stars vary either regularly or irregularly, although good statistics for the 
evolved stars likely to be SN progenitors seems to be lacking. The nature of the variability is closely related 
to the stellar type (e.g., ISzczygiel et al.ll201(J) . but not as yet in a manner that is a better diagnostic of 
the star than its location in a color magnitude diagram except in the limit of helioseismology. There are 
theoretic al arguments that the pulsational properties of red supergiants change with the onset of carbon 
burning (IHeger et al l 19971), but this has not been developed to the point of providing any observational 
guidance. Similarly, lAmett & MeakinI (1201 lb show that in three dimensions the nuclear burning fronts of 
stars in these late phases can be very dynamic, which could drive surface effects. Thus, it is likely that the 
progenitors are variable, particularly the red supergiants, but it is unknown whether the variability properties 
of SN progenitors show any recognizable difference from stars that have not commenced carbon burning. 

Another s ource of variability co mes from shells of material ejected during ou tbursts either as observed 
for SN 2006jc dPastorello et al.ll2007l) or inferred from the post-SN evolution (e.g.. lFox et"aDl20Ill) . Ejected 
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material has observable effects on the Ught curve of the progenitor if dust forms in the ejected material. 
As a dusty shell expands, its optical depth drops as r oc 1/r^ oc 1/f^ and th e star becomes steadily brighte r 
and bluer, as is observed for sources such as r] Carinae (see the review by [Humphreys & Davidsonlll994t) . 
Observing these changes, as well as the associated mid-IR emission, constrains the time of the eruption if 
unobserved and the amount of material ejected. 

Finally, close binary stars can show ellipsoidal variations from a relatively broad range of viewing 
angles and eclipses from a narrow range. For example, the progenitor of SN 1993 J had a roughly 15% 
probab ility of producing visible eclipses based on the binary evolution models of IStancliffe & Eldridge 
(l2009h . It would, however, require extraordinarily good luck to find a supernova progenitor in an eclipsing 
binary since a more typical probability is 5% for a 100% binary fraction, and for the distance limits we 
are considering (< 10 Mpc) the SN rate is only ~ 1/year . In most cases it will likely be easier to identify 
candidate binary companions once the SN has faded (see iKochanekll2009h . 



With wide field cameras on 8.5m class telescopes it is now possible to begin exploring these problems 
in nearby (< 10 Mpc) galaxies. While crowding means that the Hubble Space Telescope (HST) is generally 
needed to provide photometry, difference imaging methods make it relatively easy to monitor individual 
stars from t he ground because luminous variable stars are relatively rare and hence not crowded. This is 



illustrated in Gerke et al 



Cameras (LBC, iGiallongo et al. 



(1201 ih. where we identified over 100 Cepheids in M 81 using the Large Binocular 



2008) on the twin 8.5m Large Binocular Telescope (LBT). The shortest 



eriod Cepheids in iGerke et al.l (|2011h . at P = 10 days, have masses of order M ~ 6Mq (e.g., iBono et al. 
2000.) that are well below the mass limit for SN progenitors. We are presently monitoring 25 nearby galaxies 
in the UBVR bands using the LBT/LBC, sparsely monitoring most and more intensively monitoring a few. 
In addition to the search for variable stars, the data also enables two more speculative projects. The first 
is to set limits on the existe nce and rate of failed s upemovae, massive stars forming black holes without a 
dramatic external signature (IKochanek et all 120081) . Based on the best current statistics for s tar formation 
and supernova rates, these could represent up to half of all stellar deaths (IHoriuchi et al.ll201 ih . 



The second speculative goal is to study the varia bility of SN progenitors. Unfortunately, the first SN in 
our sample, SN 2009hd in NGC 3627 (|Monardll2009l ). occurred when we had almost n o data and it also lay 
behind a dust lane in the wings of an unobscured bright star (see lElias-Rosa et al We obtained no in- 



teresting limits on the variability of the progenitor, although we should be able to obtain UBVR photometry 
of the progenitor once the SN has faded. The second SN in our sample is SN 2011dh i n M 51 (IGriga et al. 



2011 


). A candidate r 
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The object is relatively yellow (T ~ 6000 K) and there are suggestions tha t the 



SEP may repre s ent a composite of two, presumably binary, stars. Furthermore, lArcavi et all (1201 II) and 



Soderberg et al.l (|201 ll) argue that the star which exploded must be a more compact, blue star based on the 
rapid evolution of the early-time light curve a nd spectroscopy This would be consis tent with the pres ence 
of Ha emission in pre-explosion HST images dSzczygiel et al.ll201 ih . As discussed in lKochanekl (120091) . we 
expect 50-80% of SNe to occur in stellar binaries, and it is relatively easy for the cooler star to dominate 
the optical emission. Combined with the Type lib spectroscopic type, the system seems very similar to the 
binary progenitor of SN 1993J. While we have yet to carry out the intensive monitoring phase for periodic 
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variables in M 51, we have 5/4 epochs of UBV/R data spread over three years and came close to observing 
the star just before it exploded since the next run started on 2 June 201 1 with M 5 1 as a high priority target. 
The candidate progenitor is well detected in all epochs and variable. We will simply refer to this star as the 
progenitor to avoid constant use of the clumsy phrase "candidate progenitor," but the ultimate interpretation 
of its variability clearly depends on resolving this ambiguity. In ^we describe the data and our variability 
analysis based on difference imaging. I n i3]we discuss some of the implications. We adop t a distance to 



M 51 of 8.3 Mpc (|Poznanski et al.ll2009l) and a foreground Galactic extinction of 0.035 mag (|Schlegel et al 
1998L . 



2. Data and Results 

We observed M 51 with the LBT/LBC before the SN explosion on 9 March 2008, 28 January 2009, 
19 March 2010, 11 February 2011 and 29 April 2011, and after the explosion on 5 and 9 June 2011. The 
LBC-Red camera was not available for the 11 February 2011 epoch. Table [T] summarizes the observations. 
The LBC cameras have a pixel scale of 0"224, and the galaxy was placed on the central chip #2 of each 
camera. The images were bias-corrected and fiat fielded using sky fiats following standard procedures using 
IRAF mscred tasks. Fig. [T] shows the UBVR reference images as compared to the archival HST images of 
the region. The star is clearly visible, but blended with the blue star O'.'S to the North-East. The progenitor 
dominates the BVR fluxes (based on the HST images, it represents 65%, 76% and 87% of the B, V and I 
fluxes, respectively), while the U band flux is dominated by the blue star and we see a corresponding shift 
in the location of the peak. Since the SN progenitor is not visible in the U-band images, we will not discuss 
them in great detail. 

We analyzed the images using the ISIS difference imaging package ( Alard & Lupt onlll998 . AlardlboOO ) 



with a modified star matching procedure based on Sextractor ( Bertin & AmoutsI 19961) that perf( 



lorms more 

reliably in our fields. The earlier LBC data had some significant rotations and translations relative to the 
nominal field centers, and we discovered that the standard ISIS interpolation routine mis-handles interpo- 
lation in this circumstance. Once identified, the problem was corrected by modifying the ISIS spline2.c 
routine to correctly identify the target column associated with each row. We interpolated the individual 
sub-images to an R-band reference frame using a second order polynomial for the coordinate transforma- 
tion constructed from 100-300 stars with typical rms residuals in the stellar matches of less than 0.1 pixels. 
A reference image was constructed for each band from 6-10 of the sub-images obtained prior to the SN. 
We carried out the difference imaging both on the individual sub-frames and by combining the interpolated 
sub-frames and then difference imaging the combined images. 

Fig.|2]shows the light curves with the magnitude at the epoch of highest brightness normalized to zero. 
In R band, the flux is dominated by the progenitor, and over the last three years it faded by roughly 0. 13 mag 
at a rate of 0.039 ± 0.006 mag/year. Note that for this light curve we constructed several stacked images 
for each epoch rather than combining all the data into a single image. The results for the separate images at 
each epoch are mutually consistent. The B and V light curves show an initial rise and then a fall, although 
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an increasing fraction of the light comes from the contaminating blue star at shorter wavelengths (24% and 
V and 35% at B). The U band light curve shows no convincing variability and is primarily emission from 
the contaminating star. 



We used DAOPHOT ([Stetson I Il987h to catalog the stars on ou r U, B, V and R-band ISIS reference 



images. We calibrated the U-band data using 6 standard stars from iPastorello et al.l (|2009l) and obtaining 
a calibration uncertainty 0.041 mag. The B, V, and R-band data were calibrated using a limited number 
(1 4) of Sloan Digital Sky Survey (SDSS DR-fl) stars found on chip #2 (unfortunately, the standards used 



by IPastorello et all (120091) are all satur ated in the LBT B VR images). We transformed their magnitudes 
from the ugriz to UBVRc system using iJordi et al. I (l2006l) . After eliminating stars with uncertainties larger 



than 0.3 mag, we were left with only 7, 5 and 8 stars in B, V and R, respectively, leading to calibration 
uncertainties of 0.230, 0.095 and 0.081 mag. These absolute calibrations are relatively unimportant for our 
present discussion. We then constructed light curves from the differenced images for all the sources in the 
LBT DAOPHOT catalogs. Fig.[3]shows the variance aR of these 27,000 R-band light curves as a function 
of the R-band magnitude. The overall trend with magnitude simply represents photon-counting statistics and 
we have made no major effort to clean the variability catalog of artifacts. Since ISIS tends to underestimate 
the uncertainties in light curves, we used Fig. [3] to rescale its error estimates. Given that most sources are 
not variable, the light curve uncertainties should match the median variance seen in Fig.[3](heavy line), and 
this requires scaling the ISIS light curve uncertainties upwards by a factor of three. We use these rescaled 
uncertainties in Fig. |2l The final light curve of the SN progenitor is presented in Table |2] The individual 
uncertainties are correct for the light curve, but there are additional global uncertainties from setting the 
magnitude of the source in the reference image and the zeropoints that are reported in the caption. 

Clearly we cannot interpret such a sparse light curve physically, but we can for the first time examine 
the variability of an SN progenitor as compared to other stars. Keep in mind, however, the debate on 
the identification of the progenitor - this may be the variability of a binary companion to the actual SN 
progenitor. As shown in Fig. |3l the progenitor is significantly more variable than the typical star of its 
magnitude. Since its light curve is roughly a linear decline, we made linear fits m, = st + niQ to all the light 
curves to give a slope s (mag/year, so positive slope is fading) and the residual dispersion a about the linear 
fit. In Fig.m we show the distribution of the 3800 stars within ±1 R mag of the progenitor in this space of 
ordered (s) and disordered (cr) variability. In this variability space, the progenitor properties are clearly quite 
different from the typical star, with 93% of the stars having smaller slopes in absolute value. However, note 
that simply comparing to stars of similar R band magnitude averages over stars of many types, and there are 
probably many false outliers in the distribution because we have not inspected all the light curves of these 
objects for artifacts. 

We can construct a better comparison sample by searching for "true" analogues to the progenitor in the 
B, V, I catalogs obtained from the HST images (constructed with DAOPHOT). We matched the magnitudes 
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niaj of each potential analogue to those of the progenitor mpj, 

X^= Yl {mpj-maj-AM-RiAEYao^ (1) 

allowing for a difference in luminosity \AM\ < 0.5 mag and extinction \AE\ < 0.1 and accepting those 
with < 4 for a fixed ctq = 0. 1 mag. This identified 235 such stars, of which 77 lay outside the masked 
regions of the R band LBT image. Many of the potential analogues lie in the central regions of M 51, 
which are saturated in our R band images, while the remainder trace the spiral arms. We visually examined 
the analogues to detect any potential artifacts. FigH] also shows the variability properties of these analogue 
stars (black crosses), and the slope of the progenitor is still unusually large, with only 5% (4 objects) of the 
analogs showing larger absolute slopes. 



3. Discussion 



Our primary result is that with only 4/5 epochs of ground based data there is no difficulty in detecting 
low level variability in a supernova progenitor candidate at a distance of almost 10 Mpc. In the R-band we 
find a relatively steady decline of 5^ = 0.039 ± 0.006 mag/year over a three year baseline with rms residuals 
of only 0.02 mag. With so few epochs, we cannot interpret this physically and th e uncertain t y abou t the true 
identity of the progenitor complicates the interpretation in any case. As noted by lKochanekl (|2009|), 50-80% 
of ccSNe should be in stellar binaries at the time of explosion. Furthermore, if the binary consists of a 
blue, hot star and a red (in this case yellow) cool star and the blue star explodes, it will not be uncommon 
for the visual emission from the progenitor to be dominated by the red companion rather than the star 
which exploded. This issue will be resolved as the direct emission from the SN fades, although we should 
note that the system most advanced as a possible analogue, SN 1993J, produces so much emission from 
the expanding shock moving throu gh the circum stellar medium that the binary companion only became 
observable a decade later (see Maund et al.ll2004l) . 



While the utility of finding that the progenitor is an eclipsing binary is obvious, one could legitimately 
ask whether detecting other sources of variability has any use. At a very basic level, these stars are different 
from all other stars in their galaxies. In a Hertzsprung-Russel diagram this is not apparent because there is 
nothing special about their luminosities and surface temperatures. Variability opens a new window to search 
for these differences. There may be none, but the default expectation that there are no precursor signals to 
SNe is essentially based on assuming that the only effect is the quasi-static evolution of the stars and that 
this is too slow to detect. 

How slow is the quasi-static evolution? Fig. |5] shows esti mates of the r a te of e volution for t he las t 
cent ury before collapse der ived from the models tabulated by ISchaller et al.l (Il992h . iHeger et al.l (l2000r) 
and iHeger & Langer (l2000b . The stars are generally becoming slightly fainter and hotter at rates of 0.1- 
1 millimag/year. Bear in mind that these models were not intended for this purpose. Based on photon 
statistics, it is possible to detect such slow rates of change. Given A'^ observations with an 8.5 m telescope, 
uniformly spaced over time T, with exposure time t^xp for each epoch, the 3a sensitivity to a temporal 
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gradient is 

^ 0.040 mag / D \ / (l+g)ksec \ ^^_o.2(m.+8) 

for a star of absolute luminosity My at distance D with the background flux in the extraction aperture being 
B times that of the star. For the R-band observations of M 5 1, we have T = 3. 1 years, N = 4, Ntgxp = 8000 sec, 
D = 8.3 Mpc, and B ~ 10, yielding ^'3 ~ 14 millimag/year for the My — -8 mag progenitor sta r. This result 
does not depend signific antly on the adopted SN distance - if we take D = 7.1 Mpc (used by 



Maund et al. 



20111) or D = 7.66 Mpc (IVan Dyk et al.ll201 11) . then 53 ~ 12- 13 millimag/year. The formal estimate from 
our progenitor light curve is that we reached 53 ~ 18 millimag/year, roughly consistent with this estimate 
but well above even the highest rates shown in Fig. |5] 



^Consider, however, the data we have accumulated on M 81 to study its variable stars (IGerke et al. 

20111) . With D = 3.6 Mpc, = 50, Ntej^p = 21540 sec, and T = 4.1 years, the nominal 3a sensitivity of 



^3 = 3 millimag/year (at My = -8) is below some of the quasi-static evolution rate predictions. Extending 
the time baseline to 10 years reaches 5'3 ~ 1 millimag/year, and campaigns with N = 1000 and T = 30 years 
yielding S3 ~ 0.1 millimag/year are conceivable. Thus, it is possible to obtain the photon statistics needed 
to probe this phase of stellar evolutio n, and we k now from ground-based milli-magnitude photometry of 



planet ary transits (see the review by IWinnI 120101) and ~ 10~^ milli-magnitude photometry of transits in 
space (IBorucki et al.ll2009h that systematic errors can be controlled well enough to approach the statistical 
limits. A milli-magnitude is not what it used to be. A few other regimes of stellar evolution, such as the 
post-helium flash evolution of stars onto the horizontal branch (Bildsten, private comm.) and some massive 
stars in the Hertzsprung gap should also be fast enough to observe, albeit not at Mpc distances. We do need 
theoretical models intended for making estimates of the surface evolution in these last phases, particularly 
since most studies of stellar evolution literally freeze the envelope for the last phases. 

The systematic problem that will most limit measurements of quasi-static evolution is that stars vary 
on many time scales other than this evolutionary time scale. Variability acts like an added source of noise, 
essentially adding B ~ (A/apiiot)^ to Eqn. |2j where A is the amplitude of any unmodeled variability and 
aphot is the typical photometric error at any epoch. If other sources of variability dominate, then the question 
becomes whether the variability of progenitors can be distinguished from that of other stars. With the sparse 
data we have, we can only consider simple metrics. Here we examined variability in the space of a mean 
trend and the residuals around it, finding that the progenitor has modestly unusual variability properties. 
Since there i s esse ntially no theoretical guidance on variability in these late phases, other than the study by 



Heger et al.l (Il997h that the oscillation properties change, and we are faced with the additional uncertainties 
about the nature of the progenitor, it is premature to draw conclusions. We will probably need variability 
statistics on several progenitors to begin having a clear path for interpreting the results. 

The particular properties of one ambiguous object are not a revolution. But the ability to make the 
measurements may be revolutionary. Using difference imaging techniques we can measure the variability of 
any SN progenitor within 10 Mpc from a ground-based 8.5m telescope at levels that certainly approach and 
may reach the variability expected from the quasi-static evolution of the stars. Any variability significantly 
above that level is trivial, and we simply face the quantitative question of whether the variability of post- 
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Carbon burning stars can be distinguished from that of other stars. In other words, we could ask whether 
the variability properties can be used to point at stars about to explode. The "about to" is likely on the order 
of the present duration of human civilization (the 10^ - 10^ years after Carbon ignition), but this is still a 
remarkably narrow window compared to the life times of even massive stars. The important point about our 
observations of the progenitor of SN 201 Idh is that these are now observational questions - we do not know 
if the answer will be boring or exciting, but we know we can answer the question. 



Forty galaxies produce 90% of the local (< 10 Mpc) ccSNe rate of roughly 1 SN/year (see lKochanek et al. 



2008h . It conservatively requires 4 nights on an 8m class telescope with a wide field camera (LBC on LBT, 



Suprimecam on Subaru, or to a lesser degree, IMACS/Megacam on Magellan) to provide one epoch of data 
with depth comparable to our present data for all 40 galaxies. Such single epoch data generally exists but 
is not very useful because seeing induced confusion means that photometry of individual stars is essentially 
impossible at these distances with ground-based optical resolution. Accurate single epoch fluxes require the 
high resolution of the Hubble Space Telescope (HST). Galaxies at these distances are not, however, crowded 
with l uminous variable stars even at ground-based resolutions, and, as shown in our study of Cepheids in 
M 81 (|Gerke et al.ll201 ih . variability should be measured from the ground with difference imaging and only 
the absolute calibrations done from space. Getting to the modest numbers of epochs we use here is relatively 
easy, roughly 20 nights to obtain 5 epochs for every galaxy. This is approximately where we stand in our 
LBT survey of 25 of these 40 galaxies - we have a median of 5 epochs. Reaching 30 epochs is expensive, 
but this is roughly the threshold where one can identify and phase periodic variables like Cepheids and build 
long term light curves of fainter transients. At least in our LBT survey, we are trying to reach this level 
for a subset of the galaxies that are interesting for studies of the distance scale (e.g., M 81, NGC 4258, 
M 10 1) or where there are interesting, faint transients to be monitored (e.g., SN 2008S in NGC 6946, see 



Kochanek et all (|201 ih ). Achieving the next level, 100 epochs, probably requires a dedicated imaging tele- 



scope like LSST, but only represents 5-10% of the observing time over a period of 20 years. With this many 
epochs, most ecUpsing binaries will be identified and it should also enable searches for microlensing events. 
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between the R and U-band LBT images because the U band flux is dominated by the bright blue star (easily 
seen in the HST F336W image) to the North-East. 



Table 1. Observation Log 



Date 


LBC Red 


LBC Blue 




Nr X te:cp 


Nv,NB,NuXtexp 


9 March 2008 


6x300 


2,2,2x300 


28 January 2009 


8x200 


3,2,3x200 


19 March 2010 


14 X 200 


3,6,5 x 200 


11 February 2011 




3,3,3x200 


29 April 2011 


9x200 


3,3,3 X 200 


5 June 2011 


9x200 


3,3,3x200 


9 June 2011 


9x200 


3,3,3x200 



Note. 



— Exposure times are in seconds. 
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4500 4750 5000 5250 5500 5750 

HJD - 2450000 



Fig. 2. — The U, B, V and R differential light curves of the SN 201 Idh progenitor. The mean magnitude at 
the epoch of the highest brightness was subtracted from the hght curves. There is increasing contamination 
from the nearby blue star at shorter wavelengths. The dotted line marks the date of the SN explosion. 
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18 19 20 21 22 23 24 25 26 27 

R [mag] 

Fig. 3. — Variable sources in M 51. The points show the variance cjr in the R-band light curves of ~ 27,000 
sources as a function of the R-band magnitude. The heavy line is the median of the cr« distribution within 
0.5 mag wide magnitude bins, while the two thin lines symmetrically enclose 68% of the objects around the 
median, representing an equivalent of one standard deviation in this non-gaussian distribution of ctr. The 
properties of the SN 201 Idh progenitor are marked with an open circle and we see that it lies just inside the 
upper thin hne, almost falling into the group of 16% most variable objects. We have not carefully inspected 
this sample for false sources of variability, so the significance of the progenitor variability is underestimated. 
The general trend simply represents the scahng of photon counting uncertainties with flux {gr (xR/5). 



-13- 



600 
400 
200 


0.1 

0.05 

0.01 
0.005 

0.001 
0.0005 



-0.1 -0.05 0.05 0.1 



1 — I — I — I — I — r 



T — I — I — I — I — I — m — I — I — I — I — I — I — I — I — r 



brightening 



T 

fading 



I 

I 



r „ •■■ xlO 



I I I I I! I I I I I I !l I I I I I I I 



X 




J I I I I Lj_I I I I I I L_lJ I I I I I I I I L 



J_ 



T — I — I — I — I — I — I — I — r 



random 



■1 



xlO 



_.i 



ordered 



0.1 
0.05 

0.01 
0.005 

0.001 
0.0005 



-0.1 -0.05 0.05 0.1 
slope s [mag/year] 



200 400 



Fig. 4. — Variability statistics of M 51 stars. The scatter diagram shows the distribution of objects in light 
curve slope s and rms residual a for the progenitor (large circle), for 3800 stars within 1 R mag of the 
progenitor (small points) and for 77 analogue stars of similar luminosity and spectral energy distribution 
(crosses). The dashed Unes mark the values observed for the progenitor for both signs of the slope. The 
panels to the side show projected histograms of the comparison stars, where the histogram for the small 
sample of analogues has been multiplied by ten. While the rms variability of the progenitor is typical of 
either sample, stars fading as rapidly as the progenitor are relatively rare. 
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Fig. 5. — Crude esti mates of the rate of ch ange in stellar fluxes over the last century before core collapse 
based on models from lschaller et al" ( 1992 ) (gray), Heger et al. ( 2000 ) and Heger & Langer (2000) (black), 
it None of these models were, however, intended for this purpose. A black body is assumed to convert L and 
Te into estimates of the rate of change in the B (filled square), V (open square), R (filled triangle) and I (open 
triangle) magnitudes. The different filters are slightly offset to improve visibility. The filled pentagon shows 
the slope measured for SN 201 Idh and the open pentagon below it shows our e stimated sensitivity li mit - the 
points are placed where the luminosity of the progenitor roughly matches the lSchaller et all (119921) models. 
The solid curve shows the (statistical) sensitivity of the present data at R-band, assuming N = A epochs 
spread over T = 3.1 years, with the background flux in the photometric aperture being B = 10 times that of 
the star for an exposure time t^^p = 1800 sec on an 8.5 m telescope. T he absolute V band rn agnitudes are 
estimated from lschaller et al" ( 1992 ) and the distance is set to 8.3 Mpc ( Poznanski et al. 2009). The dashed 
curve shows the sensitivity that can be achieved relatively easily from the ground, increasing the baseline to 
r = 10 years with A^^ = 10^ epochs, and the dotted curve shows a not impossible study with T = 30 years and 
N = 10^ epochs. 
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Table 2. The progenitor light curve 



Date 


HJD-2450000 


U 




B 




V 


(TV 


R 






[days] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


9 March 2008 


54534.5 


21.796 


0.050 


21.908 


0.019 


21.275 


0.024 


21.322 

21.338 


0.025 

0.027 


28 January 2009 


54859.4 


21.785 


0.050 


21.814 


0.017 


21.228 


0.027 


21.352 
21.365 


0.024 
0.031 


19 March 2010 


55274.5 


21.716 


0.042 


21.857 
21.867 


0.020 
0.017 


21.297 


0.023 


21.368 
21.369 
21.379 
21.372 


0.025 
0.025 
0.023 
0.020 


11 February 2011 


55603.5 


21.745 


0.033 


21.893 


0.012 


21.298 


0.015 






29 April 2011 


55680.4 


21.767 


0.044 


21.921 


0.012 


21.344 


0.020 


21.446 


0.018 



21.448 0.023 
21.452 0.017 



Note. — The uncertainties on the light curves do not include the uncertainties in the magnitude of the star on 
the reference image. The DAOPHOT uncertainties in the reference image magnitude are 0.033 in U, 0.051 in B, 
0.036 in V, and 0.026 in R, and the magnitude calibration uncertainties are 0.100 in U, 0.230 in B, 0.095 in V and 
0.081 in R. 
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